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ABSTRACT 

Context. It is well known that the tilt angles of active regions increase with their latitude (Joy's law). It has never been checked before, 
however, whether the average tilt angles change from one cycle to another. Flux transport models show the importance of tilt angles 
for the reversal and build up of magnetic flux at the poles which is, in turn, correlated with the strength of the next cycle. 
Aims. Here we analyse time series of tilt angle measurements and look for a possible relationship of the tilt angles with other solar 
cycle parameters, in order to glean information on the solar dynamo and to estimate their potential for predictions of solar activity. 
Methods. We employ tilt angle data from Mount Wilson and Kodaikanal observatories covering solar cycles 15 to 21. We analyse 
the latitudinal distribution of the tilt angles (Joy's law), their variation from cycle to cycle and their relationship to other solar cycle 
parameters, such as the strength (or total area covered by sunspots in a cycle), amplitude and length. 

Results. The two main results of our analysis are: 1. We find an anti-correlation between the mean normalized tilt angle of a given 
cycle and the strength (or amplitude) of that cycle, with a correlation coefficient of ?> = -0.95 (99.9% confidence level) and r^ = -0.93 
(99.76% confidence level) for Mount Wilson and Kodaikanal data, respectively. 2. The product of the cycle averaged tilt angle and 
the strength of the same cycle displays a significant correlation with the strength of the next cycle (r^ = 0.65 at 89% confidence level 
and r^ = 0.70 at 92% confidence level for Mount Wilson and Kodaikanal data, respectively). An even better correlation is obtained 
between the source term of the poloidal flux in Babcock-Leighton-type dynamos (which contains the tilt angle) and the amplitude 
of the next cycle. Further results are: We confirm the linear relationship (Joy's law) between the tilt angle and latitude with slopes 
of 0.26 and 0.28 for Mount Wilson and Kodaikanal data, respectively. In addition, we obtain good positive correlations between the 
normalized area weighted tilt angle and the length of the following cycle, whereas the strength or the amplitude of the next cycle do 
not appear to be correlated to the tilt angles of the current cycle alone. 

Conclusions. The results of this study indicate that in combination with the cycle strength, the active region tilt angles play an 
important role in building up the polar fields at cycle minimum. 
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lO ■ 1. Introduction to the poles, finally causing the reversal. It was already proposed 

^^ byl Leightonl (11969 ) that for the reversal to occur there must be 

O Solar cycles diff^er from each other, showing different lengths, cancellation between the leading portions of spots on opposite 

T-H amplitudes and strengths. Understanding the cause of such vari- hemispheres through the slight tilt of the bipolar regions. In this 

^ ■ ations and, ideally, reproducing them is one of the aims of dy- ^^y ^ ^^^^^^^ ^^^^^^^^ ^f ^^e following polarity flux reaches the 

■^ namo theory. poigg 
^^ Magnetic flux transport dynamo models of the Sun's global 

H I magnetic field have been shown to reproduce fairly well the ISchriiver et all (l2002h tested the hypothesis where the polar 
. . . . amplitude and duration, among other characteristics, of the so- magnetic field on the Sun is determined by the accumulation of 
lar cycle for a t least a few cycles (Charbonneau, 2005, 2007; field transported poleward from sunspots at lower latitudes as a 
iDikpati & Gilmanl 1200 6). Some of the key ingredients of such consequence of the tilt in the bipoles. Their model was not able 
models include differential rotation, meridional flow, latitude to reproduce the polar field measurements of the past years if 
distribution of sunspots, latitudinal drift and a systematic tilt an- only the rate at which suns pots emerge is var ied from one cy- 
gle of the bipolar groups (Joy's law). These ingredients together cle to another. Furthermore.l Wang et alJ (l2002h included a cycle 
explain the polarity reversal of the magnetic field at the poles to cycle variable meridional flow in order to achieve agreement 
every ~ 1 1 years. Due to differential rotation the magnetic field with other observations, demostrating that this variable merid- 
lines are wound up around the Sun's rotation axis and when this ional flow could serve as a regulator of the polarity reversal pro- 
field is strong eno ugh it becomes buoyant and rises to th e sur- cess. Flux-transport simulations also showed that the strength of 
face as sunspots (Babc ockLll96ltlDikpati & GiImanLl2008h . The the polar fields, which feed the dynamo and help determine the 
magnetic flux from the sunspots is carried by the meridional flow strength of the next cycle, is sensitive to the average tilt angle o f 

the active regions of the previous cycle dBaumann et aI.L |2004|) . 

Send offprint requests to: M. Dasi-Espuig Here we investigate whether there is a variation of the cycle av- 
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eraged tilt angle and Joy's law from cycle to cycle and whether 
there is a relationship between the tilt angles and the strength, 
i.e. the activity level, of the following cycle. 

Previous studies of the sunspot tilt angles have mainly 
focused on their relationship with other spot parameters such 
as magnetic flux, drift motions, rotation, area, polarity sepa- 
ration and cycle phase among others (Howard 1 199 iL Il996t 
ISivaraman et al.L 12007) . Variations from one activity cycle to 
the next, however, have never been explored. This could shed 
some light on the mechanism by which the magnetic field of 
active regions is transported to the poles and may thus have the 
potential for forecasting future solar activity. 

Prediction of future solar activity is one of the main chal- 
lenges in solar physics and is not only of scientific impor- 
tance but potentially helps to make predictions about changes 
to our natural environment that can affect our lives, e.g. space 
weather and Earth's climate. Most of the present day predic- 
tions are_based_onstatistical_anal}^se^ activity in the 
past (iHathawav et al.Lll999HHafliawavll2009l) . A more physics- 
based approach is ofifered by models of the evolution of the 
Sun's magnetic field, although recent dyna mo computations 
have given controversia l resul ts for cycle 24 (|Dik:pati & Gilmanl 
I2006t IChoudhuri et al.L l2007l iJiang et al.L I2OO7I) ! 

The paper is structured as follows. In Sect. 2 we describe the 
data, the method and some tests. Section 3 presents the results, 
which are then discussed in Sect. 4. In Sect. 5 we present our 
main conclusions. 



2. Data and Tests 

2.1. Data 

For our analysis we employ sunspot data derived from white 
light images taken at Mount Wilson and Kodaikanal observato- 
ries. These observatories have regularly observed the solar disc 
in white light since the beginning of the 20th century. The data 
we use cover the years 1917 to 1985 and 1906 to 1987 for Mount 
Wilson and Kodaikanal, respectively. This means that cycles 15 
to 21 are completely covered by the Kodaikanal (hereafter KK) 
record, but the first 4 years of cycle 15 and the last year of cy- 
cle 21 are missing in the Mount Wilson (hereafter MW) series. 
iHoward et al.l (1198 4') measured the positions and areas of indi- 
vidual sunspots on digitized MW images and then grouped the 
sunspots using a technique based on proximity. The gr ouping 
of individual sunspots was done bv' Howard et al.l (1 1 9841) by ap- 
plying a running box, 3° wide in latitude and 5° wide in longi- 
tude, centered at each spot on the solar disc. Any other spot that 
fell inside the box was included as part of the group. To distin- 
guish between the leading and following portions of the sunspot 
groups, they first computed the mass center. The portion to the 
east of the mass center was defined as the leading portion and the 
portion to the west as the following. This definition was applied 
to all sunspot g roups since they had no magnetic information 
(lHowardLll99lh . The tilt angle of a sunspot group, a, is defined 
as tan a - A0/[A/ cos 0], where 4> is the latitude of the center of 
the sunspot group and A0 and A/ are the differences in latitude 
and central meridian position between the center of gravity of 
the leading and following portions of the group, respectively. 

The final data set includes dates of observations, positions, 
area and number of individual sunspots for each sunspot group 
and for its leading and following portions, as well as the tilt an- 
gle. A description of the calculation of the tilt angles can be 
found in the paper bv iHowardI (Il99lh . The white light images 



from Kodaikanal observa tory were treated using the same tech- 
niques and procedures bv'Sivara man et all (Il993h . 

In the present study we al s o use the sunspot area data set 
compiled by Bal maceda et al.l (l2009l) using observations from 
a number of different observatories that were carefully cross- 
calibrated in order to reduce, as much as possible, systematic 
and other differences between observations at different sites as 
well as the number of data gaps. The combined data set goes 
back to 1874. 
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Fig. 1. Area versus polarity separation between the leading and 
following portions of a sunspot group for (a) Mount Wilson and 
(b) Kodaikanal data sets. The dashed line corresponds to 16°. 



2.2. Data evaluation 

The first inspection of the MW and KK data revealed a great 
number of zero values for the tilt angles (~ 22% in MW and 
~ 30% in KK). Therefore, we checked whether all zeros are 
real or just missing values. We assumed that the tilt angle can 
only be determined if the leading and following portions of a 
group have at least one spot each. Thus we neglected the zero 
tilt angle values in all cases when either portion of a group con- 
tained no spots. By applying this criterion, we found only one 
real measurement of a tilt of exactly zero degrees in the MW 
data set and none in the KK data set. The rest of the zero values 
only mark that it was not possible to measure the tilt angles for 
some reason. In order to accept a tabulated tilt angle as valid, we 
also required that the distance between the leading and following 
portions is less than 16°. This is justified by the distribution of 
sunspot group areas with polarity separation (see Fig. 1). Most 
of the groups lie in a range between - 350 microhemispheres 
in area and 0° - 16° in polarity separation. Only around 0.6% 
and 0.4% of all the groups in the MW and KK data sets, respec- 
tively, present a polarity separation bigger than 16°. Of these, 
most have areas below 70 microhemispheres and are thus most 
likely typos since their total area is small and at such polarity 
separations it is most improbable that the two polarities belong 
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to the same group. All together, i.e. based on both criteria, 
have rejected 22.5% of MW groups and 30.6% of KK ones. 
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KK data sets with a mo re complete daily sunspot data set from 
iBalmaceda et al.l (l2009h . In this data set, where it is possible to 
distinguish between gaps and spotless days, only 6.5% of the 
data are missing due to gaps in the considered period. After the 
comparison we retrieved a ~ 59% and a ~ 56% of truly missing 
dates in MW and KK data sets respectively. These remaining 
gaps seem on average to be more or less randomly distributed 
over cycle phases. They should not significantly affect the aver- 
ages over a cycle. However the gaps do affect the calculation of 
the cycle length and strength. In our work we use the length as 
defined and calculated by the National Geophysical Data Center 
(see Section 3.2). In the case of the strength, we used the daily 
sunspot area data set from Balma ceda et al.l (l2009h due to its low 
percent of gaps. To avoid systematic errors, we linearly interpo- 
lated across the sunspot area data gaps. 

Next we consider another possible source of bias. Since dur- 
ing a strong cycle more spots are found on the Sun's surface the 
grouping criterion by Howard (3° wide in latitude and 5° wide in 
longitude) could lead to an erroneous grouping for such cycles, 
because spots that do not belong to the group, but rather to a 
nearby neighbouring region, could be included. We expect such 
misclassification to occur mainly in the longitudinal direction 
due to the asymmetry of the box. It would lead to an enhance- 
ment of the longitudinal separation between the following and 
leading portions of the sunspot groups (A/ cos (p) for the strongest 
cycles. Such a spuriously increased longitude separation would 
lead to a reduced average tilt. Figure 2(a) shows the time se- 
ries of monthly means of sunspot group tilt angles weighted by 
their corresponding sunspot group areas (see right-hand Y-axis) 
for both MW (solid line) and KK (dot-dashed line) records. A 
smoothing of 4 years was necessary due to the high noise the 
data presented. 

The mean area-weighted tilt angles are calculated as follows: 



Fig. 2. In all panels: Dashed thin line and left- hand Y-axes show 
the m onthly means of the sunspot area from Bal maceda et alj 
(12009) vs. time; the solid and the dot-dashed thick lines represent 
Mount Wilson and Kodaikanal data, respectively, (a) Monthly 
area- weighted means of sunspot tilt angle smoothed over 4 years 
(right-hand Y-axis), (b) area-weighted latitude separation (right- 
hand Y-axis) and (c) area-weighted longitude separation of lead- 
ing and following parts of sunspot groups (right-hand Y-axis). 



In our analysis we considered mainly cycle averages to study 
the variaton of the tilt angles from cycle to cycle. The values of 
these averages are not greatly affected by gaps in the measure- 
ments if these are distributed homogeneously throughout each 
cycle. If, however, gaps are, for example, dominantly found in 
the ascending phase of a solar cycle the mean value of the tilt 
angle for that cycle will be lower than the real value. This is 
a direct consequence of the butterfly diagram (spots at the be- 
ginning of a cycle appear at higer latitudes) and Joy's law (tilt 
angles are higher for sunspot groups located at higher latitudes). 
Unfortunately it is not possible to distinguish between spotless 
days and gaps in the data sets of MW and KK observatories. We 
find that no spots were present or no observations were made 
on ~ 60% and ~ 55% of all days in the MW and KK data 
sets respectively. In fact most of them are located within so- 
lar activity minima, which is reasonable because the number of 
truly spotless days is largest at minimum activity. It was possi- 
ble to determine the spotless days by comparing the MW and 
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where Aj and aj are the area and the tilt angle of the sunspot 
group j, respectively. In the case of monthly means, the sum 
goes over all sunspot groups in one month, while in the case 
of cycle means, the sum goes over all sunspot groups present 
during one cycle. Also plotted are monthly means of sunspot 
area (dashed line and left-hand Y-axes) from iBaknaceda et al.l 
(2009). Cycle 19 is clearly the strongest and has, at the same 
time, the lowest values of the tilt angles of the 7 analysed cy- 
cles. We test in Fig. 2 (b) and (c) the possibility of a systematic 
error. Here we plot monthly means of sunspot latitude (A0) and 
longitude (A/ cos (f>) separations between the leading and follow- 
ing portions of sunspot groups. The solid line is again used for 
the data from MW observatory and the dot-dashed line for KK 
observatory. The drop in A0 during cycle 19 indicates that the 
low values of the tilt angles during this period are due to a lower 
latitudinal separation of both polarities while the fact that cycle 
19 is not conspicuous in Fig. 2 (c) indicates that its low tilt an- 
gles are not due to a larger longitudinal separation. This result 
does not exclude that the grouping algorithm might have com- 
bined together sunspots that with magnetic information would 
have been grouped differently. In any case we believe that, if 
there is a systematic error in the tilt angles, this is not seen in 
the latitude and longitude positions of the leading and following 
portions of sunspot groups. 
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2.3. Joy's Law 

The tilt angl e dependence on the latitude was first found by 
Joy in 1919 (Hale et al.,'1919') and later confirmed by ot her au- 
thors (iHoward. 1991; Wang & Sheelev, 1991; Sivaraman_eLalJ, 
Il999l) . It provides strong constraints on the magnetic field 
strength of the flux tubes in the tachocline, which eme rge to 
form the observed ac tive regions (D'Silva & Choudhuril Il993t 
ISchiissler et al.L [l994l) . This relation shows a positive trend: the 
tilt angles are larger for sunspots at higher latitudes. We have 
used data from MW and KK observatories in order to rederive 
this relationship as a test. In Fig. 3, the tilt angles averaged over 
the complete data sets for latitude bins of 5° are plotted versus 
latitude. MW data are represented by asterisks connected by the 
dashed line and KK observations by open circles connected by 
the solid line. The results for both data sets are in agreement 
with each other within 1 - 2c r and are very close to the results 
obtained by Sivarama n et al.l ([1999), with the mean difference 
between our points and theirs being |A| ~ 0.1° and |A| ~ 0.2° for 
MW and KK, respectively. These differences are most probably 
due to different selection criteria applied to the data. 
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Fig. 3. Mean tilt angle for bins of 5 degrees latitude vs. latitude 
for MW (asterisks connected by the dashed line) and KK (open 
circles connected by the solid line). The error bars represent ±1 
standard error. The thick dashed and solid lines are linear fits 
forced to pass through the origin to the MW and KK data, re- 
spectively. The dotted and dot-dashed lines correspond to lines 
where the slope has been taken as the ratio of tilt with latitude 
(see description in the text) for the MW and KK data respec- 
tively. 



We have also obtained linear fits to the data points excluding 
the last two bins since the number of groups in these is very low 
compared to the other bins (about 2% and 0.5% of the total) and 
the errors are higher by a factor of > 2. The fits are forced to pass 
through the origin since we expect no tilt for sunspot groups at 
the equator. The values for the fits are: 



and 

a = (0.28 + 0.06)^ 

for MW and KK data, respectively, where a represents the tilt 
angle (in degrees) and A the latitude (in degrees). The correla- 
tion coefficient of the regression lines are r^ = 0.85 for the MW 
data and r^ - 0.76 for the KK data. The linear fits are shown 
in Fig. 3, where the thick dashed line is for the MW fit and the 
thick solid line for the KK fit. Fitting the data points while taking 
into account that each point has a different value for the standard 
deviation gives the same result for the slope up to the third dec- 
imal. Also plotted are the lines a = MA where M is calculated 
as M = 'EjUj/'LjAj and aj and Aj coiTespond to the tilt angle 
and latitude of sunspot group j. The dotted line represents the 
MW data with a slope of M = 0.29 and the dot-dashed line the 
KK data with a slope of M = 0.30. This is shown as compar- 
ison since it is not affected by the fact that each bin contains a 
different number of points. 

The slopes found here a re lower than those obtained by 
IWang & Sheelevid 198911 199 lb from daily magnetograms for cy- 
cle 21 (sin 7 - 0.48 cos 6 + 0.03 where y is the tilt angle and 
the colatitude). However, our values for t he slope, 0.26 and 0.28, 
are closer to the 0.15 value deduced by ISchtissler & BaumannI 
(2006) who used a flux-transport model to fit data from MW 
and Wilcox Solar obervatories of the total photospheric field . 
The difference to the results of IWang & SheelevI ( Il989l Il991h 
could have a variety of causes, such as the different types of data 
considered (spots vs. magnetograms), differences in spatial res- 
olution (cf. Howard et al., 1984; Wang & Sheelev, 1989) com- 
bined with the dependence of tilt angle o n the size of a region 
(lD'Silva& Howard 119931: iHowardL Il993h . or the fact that they 
considered a single cycle. Our results for the slope of the regres- 
sion line show considerable scatter from cycle to cycle, even to 
the extent that we do not consider the values obtained for indi- 
vidual cycles to be particularly reliable (see Sect. 3.1). 

3. Results 

3.1. Average value of tilt angles 

It was first pointed out by ' Howard (Il99lh that the average tilt 
angle of all sunspot groups during the period 1917 - 1985 de- 
duced from MW data was 4.2° + 0.2°. For the whole MW data 
set period (1917 - 1985) w e obtained a va lue of 4.25° + 0.18°, 
in excellent agreement with lHowardl(ll991h . while for the whole 
KK data set (1906 - 1987), we deduced 4.51° + 0.18°. 



Table 1. Area-weighted mean tilt angles in degrees for each 
cycle for MW and KK records. 



Cycle 


MW± Icr 


KK± iCT- 


15 


5.69 ± 0.57 


5.00 + 0.50 


16 


5.08 ± 0.46 


5.91 ± 0.43 


17 


5.83 ± 0.42 


6.41 ± 0.41 


18 


5.69 ± 0.35 


4.97 ± 0.38 


19 


3.84 ± 0.33 


4.59 ± 0.38 


20 


4.63 ± 0.38 


5.73 ± 0.36 


21 


5.30 ± 0.40 


5.37 ± 0.42 



a = (0.26 + O.Q5)A 



Next we treat MW and KK data on a cycle-by-cycle basis, 
obtaining a different value of the average sunspot tilt angle for 
each cycle. Figure 2(a) displays monthly area-weighted means 
of sunspot tilt angles smoothed over 4 years through cycles 15 
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to 21 for the MW and KK data sets. Table 1 gives area-weighted 
cycle averages and Icr standard eiTor. Note that cycles 15 and 
21 are not complete in the MW data set, as discussed in Sect. 
2.1, and thus the value of the mean tilt angle for cycle 15 might 
be underestimated and the value for cycle 21 overestimated, ac- 
cording to the combination of Maunder's butterfly diagram and 
Joy's law. The low value of the sunspot tilt angles for cycle 19, 
as compared to the other cycles, indicated by Fig. 2a, is also seen 
in the cycle averaged values. 

In addition to the average tilt angles, we also attempted to 
determine Joy's law per cycle in the same manner as done in 
Sect. 2.3. However, the scatter of the individual values of the 
mean tilt angle per 5° latitude bins turn out to be too large. The 
correlation coefficient for a linear regression to the points are for 
some cycles as low as p- - 0.17 for MW and r^ = 0.026 for 
KK. Also, the errors in the calculated slopes are comparable to 
or slightly bigger than the diff'erence between these values from 
MW to KK data sets. Therefore, no clear difference could be 
determined between the slopes of Joy's law from cycle to cycle. 



3.2. Cycle parameter definitions 



For the parameter study we focus on three main characteris- 
tics of a solar cycle: strength, amplitude and length. Strength 
is defined as the total surface area covered by sunspots through- 
out a given solar cycle. We calculate it from the daily sunspot 
area data set compiled by Balmacedaet al. (2009) as the inte- 
gral of sunspot area over the duration of each cycle. This record 
is used since it has significantly fewer data gaps than the MW 
and KK data sets, as discussed in Sect. 2.2. The cycle ampli- 
tude is the highest value of monthly averaged sunspot number 
and the length is the period of time between two consecutive 
minima. Times of solar activity minimum, amplitudes and the 
lengths of cycles are taken from the National Geophysical Data 
Center; http://www.ngdc.noaa.gov/stp/SOLAR/getdata.html 

We looked for possible relationships of these parameters 
with four different quantities based on the tilt angles: cycle mean 
tilt angle, (a), cycle mean tilt angle normalized by the mean lat- 
itude of sunspots during that cycle, {a}/{A}, cycle mean area- 
weighted tilt angle, (a^^), and the cycle mean area-weighted 
tilt angle normalized by the mean latitude of sunspots dur- 
ing the same cycle, {»„)/(/(). (For a brief discussion of how 
these choices are influenced by the scatter in the tilt angles see 
Appendix A). The area-weighted tilt angles are used to give 
more importance to the bigger groups, which exhibit less scatter, 
and the normalized tilt angles are considered in order to remove 
the effect of the latitudinal dependence (Joy's law) on the cycle- 
averaged (area-weighted) tilt angles. Note that for the MW data 
set, cycles 15 and 21 are not taken into account in the relation- 
ships concerning {a} and {a„} due to their incompleteness and 
could be thus biased by Joy's law. This is not the case for the 
quantities {a)/{A} and {a^}/{A} since normalizing by the mean 
latitude removes this sourc e of bias. Sunspots in stronger cy- 
cles lie at higher latitudes ( Solankiet al.L |2008'), so that simply 
due to Joy's law these cycles would have larger mean tilt angles. 
Dividing by the mean latitude largely removes this difference 
(both, Joy's law and the dependence of mean latitude on cycle 
strength are linear), so that {a}/ (A) and {aJ)l{A) indicate intrin- 
sic changes of Joy's law from cycle to cycle. 



3.3. Relationships witiiin the same cycle 

We first investigate the possible relationship of the cycle aver- 
aged sunspot tilt angles with the three solar cycle parameters of 
the same cycle. These relations may help to shed light on the un- 
derlying magnetic flux tubes at the base of the convection zone 
and the processes that affect them on their way to the surface 
(in the case of the strength and amplitude of the cycle) and on 
the possibility that the tilt angles of active regions are involved, 
along with other features (e.g. meridional flow), in the regula- 
tion of the cycle period of the dynamo (in the case of length), or 
conversely are influenced by it. 
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Fig. 4. Cycle averaged tilt angle normalized by the emergence 
latitude vs. strength of the same cycle. The error bars represent 
Icr errors and the dashed line is a linear fit to the points. Panel 
(a) displays the results based on MW data (r^- - -0.95) , where 
cycles 15 and 21 are shown as squares and dashed lines for the 
error bars, and panel (b) on the KK data set (r^ - -0.93). 



We calculated linear correlation coefficients between the 3 
solar cycle global parameters and the 4 quantities based on the 
tilt angles (see Sect. 3.2). Due to the low number of cycles, we 
also determined the probability that the correlations are due to 
chance (P). These are calculated from the probability density 
function of the student's t-distribution, which depends both on 
the correlation coefficient and the number of points in the sam- 
ple. All the values are listed in Table 2 for MW and KK data. 
Table 2 suggests that both the strength and the amplitude of a 
cycle show a significant negative correlation with the average 
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Table 2. Correlation coefficients between the 4 quantities based on the tilt angle and the strength (S ), amplitude (A) and length (L) 
of the same cycle. 



Mount Wilson 


Kodaikanal 


Parameter 




S 




A 


L 






5 




A 


L 






re 


P 


re 


P 


re 


P 


re 


P 


re 


P 


re 


P 


(a) 


-0.59 
-0.48 


0.30 
0.41 


-0.60 
-0.48 


0.28 
0.41 


-0.29 
-0.46 


0.64 
0.44 


-0.11 
-0.46 


0.04 
0.30 


-0.69 
-0.66 


0.09 
0.11 


-0.58 
0.19 


0.17 
0.68 


{a)l(A) 


-0.95 
-0.81 


1 X 10-^' 
0.03 


-0.83 
-0.91 


0.02 
4 X 10-3 


-0.40 
0.08 


0.37 
0.86 


-0.93 
-0.80 


2 X 10-" 
0.03 


-0.82 
-0.91 


0.02 
4 X 10-3 


-0.48 
0.03 


0.30 
0.95 



Correlation coefficients are represented by r^, and the probability that the correlation is due to chance by P for both the MW and KK data sets. 



tilt of the same cycle, {a), for at least KK data. For MW data, 
the probabilities that the correlations are due to chance, P, are 
about 30%, but for KK data (that includes both cycles 15 and 
21), the corresponding probabilities are lower than 10%. These 
correlations are significantly strengthened once we eliminate the 
enhanced effect of Joy's law on cycles with sunspots on average 
at higher latitudes by considering (a) /{A}. The probabilities then 
fall to values below 2% for both MW and KK data sets. For the 
area-weighted tilt angles, the correlation coefficients are weaker. 
Although these are also strengthened after the normalization by 
(A), reaching probability values below 3%, they remain slightly 
higher than for (a) /{A}. The correlations between the length and 
the 4 tilt angle based parameters are in general low, of low con- 
fidence and inconsistent in sign between the two data sets. 

Figure 4 shows (a',)/(/l,) versus S ,, where / is the cycle num- 
ber. The dashed line represents a linear fit to the points and the 
error bars correspond to Icr errors calculated by means of er- 
ror propagation, where the errors for the mean tilt angle and the 
mean latitude coiTespond to their standard error. The error bars 
have been calculated assuming Gaussian statistics and are thus 
overestimated. In MW data (Fig. 4(a)) cycles 15 and 21 are rep- 
resented by squares and dashed lines for the error bars to denote 
their incompletenes. Note that all data points lie roughly within 
Icr of the regression lines. This suggests that given the accuracy 
of the measured tilt angles (given largely by the scatter shown by 
active regions) the obtained coiTelation coefficients are near the 
maximum value achievable for data with such large uncertainty. 



3.4. Relationships witii the following cycle 

Prediction of future solar activity is important not only for space 
weather and climate, but also to test cuiTent dynamo models. In 
this section we investigate how the cycle averaged tilt angles are 
related to the global parameters of the next cycle. We calculated 
the correlation coefficients between the tilt-angle parameters (a), 
{a) I {A), (aoj) and (a'a,)/(/l) of cycle / and the parameters S , A 
and L of cycle i + 1 and the probability that these coiTelations are 
due to chance. 



Table 3. CoiTelation coefficients between the 4 quantities based 
on the tilt angle and the length, L, of the next cycle. 



Mount Wilson 


Kodaikanal 




re 


P 


P 


(a) 


-0.88 
-0.77 


0.05 
0.13 


-0.32 0.48 
-0.57 0.18 


{a)/{A) 
{aJ/{A) 


-0.46 
-0.67 


0.30 
0.10 


-0.37 0.41 
-0.61 0.15 



Correlation coefficients are represented by r^ and the probability that the 
correlation is due to chance by P for both the MW and KK data sets. 



In general the coiTelations of the strength and amplitude with 
the 4 tilt angle based quantities are low and inconsistent be- 
tween the two data sets. Only correlations between ((Q'n)/(/l))i 
and Li+i appear to be statistically significant for both data sets. 
Table 3 lists the correlations found between the tilt angle based 
parameters and the length of the next cycle. For MW data, the 
probabilities that the correlations are due to chance are below or 
around 10% for all the averages considered except for {a)/{A). 
The strongest correlation found is with the mean tilt angle ((a)) 
of value Tc - -0.88. In contrast, for KK data only {a^)l{A) 
presents a correlation {jc - -0.61) with a chance probability 
below 15%. This suggests that if the tilt angles are large, then 
the next cycle will be short. The fact that the correlation of 
tilt with the length of the next cycle is significant, but is poor 
with the strength of the next cycle is consistent with the find- 
ing that the le ngth and strength of a cycle are poorly correlated 
(r, ^ -0.37 ICharbonneau& Dikpatil 120001 and r. = -0.35, 
ISolanki et al.ll2002h 

Now, tilt angle s influence the amount of magnetic flux reach- 
ing the poles ( Baumann et all 1200 4) and the polar magnetic flux 
during activity minimum has been found to be o ne of the proxies 
that best predicts the strength of the next cycle dMakarov et all 
1989; Dikpati et al., 2008). However, the tilt angle is not the only 
parameter influencing the polar flux, which is in line with the 
poor correlation found between tilt angles and the strength of 
the following cycle (values range from 0.40 to 0.54 for MW 
and -0.58 to 0.19 for KK and are not even consistent in sign 
between the two data sets). Obviously, the total amount of mag- 
netic flux emerging over a cycl e, 4)tot, is another cent ral param- 
eter influencing the polar flux ('Bau mann et al.L 120041) . Hence a 
more physically motivated quantity to consider is <ptot{o:) / {A} or 
<f>tot{<Xw}/{A}. Since no regular and consistent magnetic informa- 
tion is available prior to cycle 20 we use sunspot areas as prox- 
ies of (plot- I.e. instead of (ptot{a)/{A) we determine S{a}/{A}. 
Sunspot areas are proportional to the amount of magnetic flux 
emerging through the s pots since the fie ld strength averaged over 
a sunspot is similar (ISolanki & SchmidtLll993l) . 

Figure 5 shows (5(q;h.)/(/1)),' versus 5,+i. Again the er- 
rors are treated by means of error propagation and assuming 
Gaussian statistics and are thus overestimated. In this case we 
were not able to propagate the errors precisely since we have 
no information on the errors of the individual measured sunspot 
areas that would affect the calculation of 5, or 5,+i. The error 
bars in Fig. 5 are calculated assuming that S j is known accu- 
rately. Using S instead of (f),ot also means that we implicitely 
assume that the ratio of magnetic flux in faculae and network 
to that in sunspots is the same for each cycle. Consequently, the 
plotted 1 cr error bars represent lower limits to the true uncertain- 
ties. Both data sets show a moderate positive correlation between 
(S {aii)l{A))i and S ,+i (see upper row of Table 4 for r^- and P val- 
ues). It is interesting to point out that both data sets, although in- 
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dependent, display almost identical fits: y, - 0.205, +1+495066.6 
and y, = 0.205 ,+i + 535359.8 for MW and KK data sets respec- 
tiely, where y = S {a„)/{A}. 
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Fig. 5. Strength of cycle multiplied by normalized mean area- 
weighted tilt angle vs. the strength of the next cycle for (a) MW 
data (re = 0.65) and (b) KK data (r^ = 0.70). The dashed lines 
are linear fits to the data points and the error bars represent Icr 
error (assuming Gaussian statistics). For the MW data set cycles 
15 and 21 are shown as squares and dashed lines for the error 
bars to indicate their incompleteness. 

Finally, we carry out a variant of the above analysis that 
is guided by dynamo models that include th e influence of the 
merid i onal circulation at the solar surface (IWang & Sheelevl 
Il99lt IChoudhuri et al.L Il995h . According to such models the 
amount of flux reaching the poles depends (for a fixed differ- 
ential rotation, meridional circulation and diffusion rate) on the 
tilt angles and the latitude distribution of the active regions. If a 
region is at relatively high latitude, then in general both polari- 
ties are dragged to the pole by meridional circulation, leading to 
a negligible change in the magnetic flux there. For active regions 
close to the equator the magnetic flux of the leading portion can 
reach and cancel the opposite polarity of the leading portion of 
an active region in the other hemisphere. This leads to an im- 
balance in the sense that mainly flux from the following polar- 
ity reaches the pole. Hence regions at low latitudes contribute 
disproportionately to the reversal and accumulation of magnetic 
field in the poles. This is thought to affect the strength of the next 
cycle since the polar fields are the input for the next cycle. We 



take into account this latitude dependence by multiplying an ex- 
ponential function of the latitude to the area- weighted tilt angles. 
The monthly means of area and latitude weighted tilt angles are 
computed as follows: 



aa,A - 






where Ay is the area of the sunspot group j from MW and KK 
data sets, Oj the tilt angle of the same group, Aj its latitude and Aq 
is a constant that determines how rapidly the exponential func- 
tion drops with latitude. The value of Aq depends on the lati- 
tudinal velocity profile of the meridional flow, which (for rea- 
sons of symmetry) is zero at the equator Small values of Aq 
correspond to a meridional flow whose horizontal component 
increases rapidly with A. In the absence of clear observational 
constraints we have set Aq to 10°, 20° and 30° and have repeated 
the analysis for each of these values. 

In Fig. 6 we plot Sa^ (solid curve) for the whole time series 
smoothed over 24 months, with S being t he monthly means of 
sunspot area froml Balmaceda et al.l (l2009h . The dashed curve is 
the 12 month smoothed S and the solid curve has been shifted by 
-1-11 years to better compare Sa^ with S of the following cycle. 
Since the whole curve has been shifted by a constant value and 
each cycle has a different duration, the lengths of the cycles of 
the solid and dashed curves do not match. However, it is seen that 
the consecutive rise in strength from cycles 16 to 19 and the drop 
from cycle 19 to 20 are reproduced. The correlation coefficients 
between the peaks of both curves reach a maximum value of 0.79 
(P= 3%) for MW and 0.78 (P=4%) for KK when taking ^o = 10° 
(see lower row of Table 4). All of the 7 cycles are considered in 
both data sets since the maxima of cycles 15 and 21 are included 
in the MW record. The correlation coefficient values range from 
0.74 to 0.79 for the MW data set and 0.78 to 0.79 for the KK 
data set when smoothing over 24, 36 and 48 months. Since the 
first 4 years of cycle 15 are not complete in the MW data set, 
we chose a 24 month smoothing as optimal to reduce the noise 
while not losing the maximum of cycle 15. 

Table 4. Correlation coefficients between expressions containing 
a of cycle / and the strength or maximum amplitude of cycle i+l. 



Mount Wilson 


Kodaikanal 


r, P 


fc P 


_(^5),vs.S,+ i 
maxiSaaj) + liyrs vs. max{S) 


0.65 0.11 
0.79 0.03 


0.70 0.08 
0.78 0.04 



Correlation coefficients are represented by r^ and the probability that 
the correlation is due to chance by P for both the MW and KK data 
sets. The two rows correspond to different expressions explained in the 
main text. 



4. Discussion 

Our understanding of the solar dynamo remains i ncomplete de- 
spite the large amount of effort invested into it (ICharbonneaui 
2005.) . One hindrance to a better understanding is the limited 
number o f observationa l constraints, some of w hich are re- 
viewed bv lGilmanI (119861 I2OO21) and Rempel|(|2008|). 

In this paper we have analyzed the time-dependence of the 
tilt angle, which has added two more observational constraints 
that dynamo models must satisfy. The first is that there is an 
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Fig. 6. Comparison of actual and predicted sunspot area. The 
dashed curv e shows monthly sunspot areas smoothed over 12 
months from Balmaceda et alJ (12009.) . The solid curve is the pre- 
diction based on the tilt angles and sunspot areas of the previous 
cycle for (a) MW and (b) KK records, both smoothed over 24 
months. 



inverse correlation between the strength of a cycle and the tilt 
angle (a) of sunspot groups observed during that cycle. This cor- 
relation was found to increase when the latitudinal dependence 
of Joy's law was taken into account (that is when alA was con- 
sidered instead of a). 

The results in the previous section are based on observations. 
While this paper concentrates on the observational signature, it 
is worthwhile to digress and speculate on the possible causes of 
this inverse correlation. In doing so we emphasise that our dis- 
cussion is only speculative whereas the results of the previous 
section are based firmly in the observational data sets. One possi- 
bility is that the field strength of magnetic flux tubes in the over- 
shoot region below the convection zone is larger during strong 
cycles. Since stronger flux tubes are less affected by the Coriolis 
force this would explain the observed correlation. In Babcock- 
Leighton type dynamos the toroidal flux tubes at the base of the 
convection zone are believed to be the result of the differential 
rotation winding up the poloidal magnetic field. The magnetic 
energy density of the loops formed in this way is likely to be 
limited to equipartition values with the kinetic energy density 
of the differential rotation. This gives a magnetic field strength 
of s» 10"^ G. Such a loop can then loose mass via an instabil- 
ity which drains mass from the slightly sub-adiabatic region 
where the flux tube is located into the covection zone, which 
increases the field stre ngth to «i 10^ G over a timespan of ap- 
proximately 6 months ( iRempel & Schusslerll2001.) . As the flux 
tube becomes stronger, it becomes subject to the Parker instabil- 
ity which causes it to erupt to the surface. The Parker instability 
depends on both the field strength and on the sub-adiabaticity of 
the layer where the flux tubes are stored: magnetic fields of tubes 
that are stored slightly deeper can become stronger prior to the 
onset of the Parker instability. 

To explain the observations requires either that in strong 
cycles the tubes ar e produced or stored slightly deeper 
dCalipari et al.L Il995h . or the region where the flux tubes are 
stored becomes slightly more stable, or that the intensification 



process acts more quickly so that higher field strengths can be 
reached before the tubes erupt. A combination of these processes 
may also be at work. 

There are a number of nonlinear, competing, factors which 
are likely to be relevant. For example the increased flux of a 
strong cycle will be more resistent to downward pumping, and 
will perhaps decrease the depth of the convective overshooting at 
the base of the convection zone. This might decrease the depth 
at which the flux is located, but will also affect the thermody- 
namic properties of the layer. Another effect, which acts in the 
correct direction, is the magnetic tension associated with the en- 
hanced poloidal flux of strong cycles. This will tend to pull the 
field lines deeper into the overshoot region, however the effect is 
likely to be weak, possibly depending on how the poloidal field 
is structured. 

In relation to changing the subadiabaticity, 

Rem pel & Schiisslen (l200lh have argued that the energy to 
intensify the toroidal flux tubes to lO^G comes from moving 
material along the entropy gradient near where the tubes are 
stored. The amount of energy involved has been estimated (e.g. 
Steiner, 20 04) to be approx i matel y lO"*" ergs. On the observa- 
tional side. iBaldner & Basul (l2008l) reported a 10"^% change in 
the wave speed squared near the base of the convection zone. 
The observed change is small and anticorrelated with activity. 
Its effect in the current context is to change the subadiabaticity 
- enhancing the stability in the region where the flux tubes are 
stored. How strong the effect would be, and how it balances 
with other effects, needs to be evaluated. 

Another possibility is that the observed tilt angles have 
been influenced by the near-surface flows. These inflows con- 
sist o f a time-dependent compon ent of the solar differential rota- 
tion dHoward & Labonta Il980h called zonal flows, and a time- 
dependent component in the meridional plane which has been 
observed by tracking magnetic featur es (Komm et al., 1993) and 
with helioseismology ( Basu & A ntia. 2003). Some models sug- 
gest that both components of th e inflow are driven by t he exc ess 
cooling associated with plage dSpruitl I2003L iRempelL l2006l) . If 
so, the strength of the inflow will be determined by the amount 
of plage, which is directly related to the strength of the current 
cycle. 

Both components of the inflows will tend to decrease the tilt 
angle over time. This effect acts on the flux tube both as it rises 
through the inflow and as it evolves on the surface after emer- 
gence. The sign of this effect is correct, enhanced inflows during 
strong cycles will reduce the observed tilt angles. To estimate 
the magnitude of the effect we now concentrate on the merid- 
ional component of the inflow (the effect of the time variations 
of the zonal flows turns out to have a similar magnitude). 

The time dependent meri dional flow includes an infl ow to- 
wards the active region belts (IZhao & Kosovichevi l2004l) which 
had an amplitude of +5m/s for cycle 23. To estimate the expected 
magnitude of the effect on the tilt angle, we assume that sunspots 
are, on average, subject to this flow for ^ 5 days before they are 
observed. These 5 days include the rise time through the flow, 
the sunspot formation time as well as a delay caused by the fact 
that not all sunspots appear on the side of the Sun facing the 
Earth. The maximum relative velocity with which the two po- 
larities could be driven towards one another by such a flow is 
lOm/s. This maximum is unlikely to be reached, so for the pur- 
poses of obtaining a preliminary estimate we assume that they 
actually move towards each other with half this speed, that is 
5m/s. Over the course of the 5 days this gives a 2. 16Mm decrease 
in the latitudinal separation of their leading and trailing fluxes. 
For an active region with a longitudinal separation between the 
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leading and trailing fluxes of lOOMm, this is a decrease in tilt 
angle of around 1.23 degrees. If we allow for the fact that cy- 
cle 19 was stronger than cycle 23, so that its inflows would have 
been stronger, then the magnitude of the effect is approximately 
consistent with the observations. 

The above are not the only possibilities for the observed neg- 
ative correlation, and have been presented only to give an indica- 
tion of some of the different types of possibilities. It is possible 
that some of these explanations will be able to be excluded on 
either observational or theoretical grounds. More modelling and 
observations will be required to pinpoint the main mechanism. 

The second result is that there is a reasonably strong cor- 
relation between the product of the strength of a cycle and its 
average tilt angle and the strength of the next cycle, r^ - 0.65 
and re = 0.70 for the two data sets, respectively. This product 
was considered because it corresponds to the poloidal source 
term in dynamos based on the Babcock-Leighton idea. The cor- 
relations were found to improve when the quantities were made 
to more closely match the poloidal source term of the models. 
Specifically, we found an improvement when we included a lat- 
itudinal dependence designed to model the effectiveness of flux 
emergence in producing global poloidal fields (which depends 
upon some of the flux crossing the equator so that regions emerg- 
ing close to the equator are more effective). This observational 
constraint supports the flux- transport dynamo model. It shows 
that the strength of a cycle is correlated at the 79% or 78% level 
(depending on the data set) with the poloidal flux of the previ- 
ous cycle. Importantly, the drop in strength from cycle 19 to 20 
is well reproduced. This tilt angle then appears to account for a 
substantial part of the variation from cycle to cycle of the activity 
level. 

This is good news in two regards. Firstly, it suggests that, 
by measuring the tilt angle and amplitude of a cycle, we will be 
able to make early predictions of the strength of the proceeding 
cycle. The predictive accuracy is not higher than, f or ins tance 
methods based on precursors (see ^ Hathawav et al.L Il999h . but 
can be made much earlier. This can be seen in Figure 6 where 
the "predictions" have been shifted by eleven years (note that 
a two year smoothing has been performed which reduces the 
predictive horizon to 10 years). A possible improvement of pre- 
dictive skill at later times might be possible by combining dif- 
ferent schemes, but this will depend on how independent they 
are. Secondly, it suggests that a major part of the nonlinear cy- 
cle modulation is associated with the tilt angle. Several possible 
non-linearities were discussed above, such as the near-surface 
inflows, the depth at which the tubes are stored and the proper- 
ties of the plasma near the base of the convection zone. More 
work is required to distinguish between these and other possibil- 
ities. 



5. Conclusions 

We have analyzed the sunspot data from Mount Wilson (MW) 
and Kodaikanal (KK) observatories in order to study Joy's law, 
the variation of sunspot group tilt angles from cycle to cycle and 
the relationship of this variation with 3 solar cycle parameters: 
strength, amplitude and length. The correlations found from the 
analysis are listed in Tables 1, 2, 3 and 4. From the analysis we 
highlight the following: 

(1) A Hnear fit to Joy's law gives a - (0.26 + 0.05)/l for the 
MW and a = (0.28 + OM)A for the KK data sets. Here a is the 
tilt angle and A the latitude, both expressed in degrees. 



(2) The mean tilt angle changes from cycle to cycle (Fig. 2a 
and Table 1). The range of values exceeds the uncertainties in 
the cycle-averaged tilt angles. 

(3) A negative correlation between the normalized tilt angle, 
or {a)l{A), and the strength of the same cycle is found (r^ = 
-0.95 and r^ = -0.93 for MW and KK data sets, respectively). 

(4) We also find a negative correlation between the latitude 
normahzed area weighted tilt angle ((a'„,)/(/l)) and the length of 
the next cycle (r^ - -0.67 and Kc - -0.61 for MW and KK data 
sets, respectively). 

(5) Finally, we discovered a positive correlation between 
the strength of one cycle multiplied by its mean area- and 
latitude- weighted tilt angle, (S {a„}/{A})i, and the strength of 
the next cycle, St+i, ijc - 0.65 and r^ - 0.70 for MW and 
KK data sets, respectively). Higher correlation coefficients are 
obtained between a tilt-angle based expression obtained through 
guidance from Babcock-Leighton type dynamo models and the 
amplitude of the next cycle (r^ = 0.79 and r^. - 0.78 for MW 
and KK, respectively). 

These results show the importance of the tilt angle of sunspot 
groups for both the prediction of solar activity and the under- 
standing of the physics behind the solar dynamo. 
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This result indicates that calculating the cycle-to-cycle devia- 
tions from a reference Joy's law is identical to calculating Joy's 
law for each cycle and then subtracting a fixed constant. A sim- 
ilar result can be shown if e, is assumed to be dependent on the 
area, as was considered in the main text. 



Appendix A: Determining the cycle-to-cycle 
variaitons in the presence of a large intrinsic 
scatter. 

As has been metioned, the tilt angles of individual active regions 
are largely random, with Joy's law being a relatively small bias. 
The purpose of this appendix is to discuss the calculation of 
cycle-to-cycle changes in Joy's law from the data. For this pur- 
pose we will assume that Joy's law applies for each cycle and 
that the scatter in the data is random and unbiased. More explic- 
itly we assume 

1. the tilt angle, a,, for each spot, /, obeys Joy's law 

a,- = a„/l; -1- 6; (A.l) 

where a„ is the (possibly cycle-to-cycle dependent) constant 
of proportionality for cycle n and e, represents the random 
deviation from Joy's law of individual sunspot groups. 

2. The e, are independent realizations of a random process with 
a mean of zero. 



We calculate our estimate, b„, ofa„ for each cycle according 



to 



b„ 






(A.2) 



where the sum is again over spots in cycle n. The error of the 
approximation is 



^n ^n ~ ^n 






(A3) 



We also consider a, the value of a based on the whole data 
set, ignoring cycle-to-cycle changes in a. The equivalent esti- 
mate, b of a has the summation extended to all cycles. Note that 
we could have also considered the cycle-to-cycle deviation of 
Joy's law, dn, from the estimate obtained over all cycles b. This 
however has exactly the same error as does «„ as can easily be 
seen: 









d„ = a-b + 

Li 'ii 



(A.4) 
(A.5) 

(A.6) 

(A.7) 



